Quantum 120° model on pyrochlore lattice: orbital ordering in MnV204 
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We present an analytical model of orbital ordering in vanadium spinel MnV2 04. The model is 
based on recent first-principles calculation indicating a strong trigonal distortion at the vanadium 
sites of this compound [Phys. Rev. Lett. 102, 216405 (2009)]. At the single-ion level, the trigonal 
crystal field leaves a doubly degenerate atomic ground state and breaks the approximate rotational 
symmetry of t2g orbitals. We find that the effective interaction between the low-energy doublets 
is described by a quantum antiferromagnetic 120° model on the pyrochlore lattice. We obtain the 
classical ground state and show its stability against quantum fluctuations. The corresponding orbital 
order consisting of two inequivalent orbital chains is consistent with the experimentally observed 
tetragonal symmetry. A periodic modulation of electron density function along orbital chains is 
shown to arise from the staggering of local trigonal axes. In the presence of orbital order, single-ion 
spin anisotropy arising from relativistic spin-orbit interaction stabilizes the experimentally observed 
orthogonal magnetic structure. 



I. INTRODUCTION 

Geometrically frustrated magnets with orbital degen- 
eracy exhibit a variety of complex ground states with un- 
usual magnetic and orbital ordersii"— Not only do these 
spin-orbital models deepen our understanding of systems 
with competing degrees of freedom, they also describe 
the low-energy physics of several transition-metal com- 
pounds. Of particular interest is spin-orbital model on 
three-dimensional pyrochlore lattice where geometri- 
cal frustration between nearest-neighbor spins leads to 
a macroscopic degeneracy in the classical ground state. ^ 
Experimentally, a signature of strong frustration is the 
occurrence of a magnetic phase transition well below the 
Curie- Weiss temperature. For systems with degenerate 
orbitals, magnetic frustration is partially relieved in the 
presence of a long-range orbital order, which is usually ac- 
companied by a simultaneous structural distortion due to 
Jahn- Teller effect. The resulting magnetic order depends 
critically on the orbital configuration and the details of 
spin-orbital interactions. 

Recently, much attention has been focused on vana- 
dium spinels AV2O4, where vanadium ions occupying the 
structural i?-sites form a pyrochlore latticei^"— The two 
d electrons of V"^"*" ion have a total spin S' = 1 and oc- 
cupy two out of three t2g orbitals. Thermodynamically, 
vanadium spinels with a non-magnetic A-site ion [A ~ 
Zn, Cd, and Mg) exhibit similar behavior: a structural 
transition which lowers the crystal symmetry from cubic 
FdSm to tetragonal lAi/amd is followed by a magnetic 
ordering at a lower temperature)^^— Assuming that xy 
orbital is occupied at all sites due to flattened VOg octa- 
hedra, minimization of a Kugel-Khomskii type Hamilto- 
nian on pyrochlore lattice gives rise to a staggered order- 
ing of the remaining yz and zx orbitals (so-called ^-type 
ordering). IS The resulting symmetry lAi/a, however, is 
incompatible with experimental observations. 

On the other hand, assuming a large relativistic spin- 
orbit coupling, a ferro-orbital order in which one electron 



occupies the low-energy xy orbital, whereas the other one 
is in states \yz)zLi\zx) has been proposed in Ref. 17'. The 
occurrence of complex orbitals imply a nonzero orbital 
angular momentum. This model successfully explains 
important experimental results: the uniform occupation 
of orbitals is compatible with space group lAi/amd; an 
ordering of orbital moment opposite to local spin direc- 
tion is also consistent with the observed reduced vana- 
dium moment. A ground state with complex orbitals 
has also been confirmed by mean-field^® and ah initio 
calculationsj ^^'^" 

Interest in antiferro-orbital order is rekindled by a 
recent experimental characterization of another vana- 
dium spinel MnV2 04r^"— where the A-site Mn^+ ion 
is in a 2>d^ high spin configuration [S = 5/2). In con- 
trast to other vanadium spinels, MnV204 first under- 
goes a magnetic transition at Tp = 56 K into a collinear 
ferromagnetic phase with Mn and V moments aligned 
antiparallel to each other. At a slightly lower temper- 
ature Tg = 53 K, a structural distortion lowering the 
crystal symmetry to tetragonal /4i/a is accompanied by 
an ordering of the transverse components of V spins. 
The ground-state orbital configuration is suggested to 
be the A-type antiferro-orbital order— Also contrary to 
collinear magnetic order in other vanadium spinels, a pe- 
culiar non-collinear order with transverse component of 
vanadium spins forming an orthogonal structure in the 
ah plane was observed in MnV2 04— 

Recently we have demonstrated the stability of orthog- 
onal magnetic structure in the limit of strong relativis- 
tic spin-orbit coupling— However, findings from first- 
principles calculation indicate a significant trigonal dis- 
tortion at the vanadium sites of MnV204,^^ whose effect 
is yet to be understood. The same authors find an orbital 
order consisting of two inequivalent orbital chains similar 
to the A-type order. More importantly, they observe an 
additional modulation of electron density profile within 
each orbital chain: the orbitals rotate alternatively by 
about 45° along the chain. This complex orbital pattern 
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is also supported by a recent NMR measurement i« 

Based on these observations, we present an analytical 
model of spinel MnV204 assuming that the t2g orbitals 
is split into a singlet and a doublet by a strong trigonal 
crystal field. As one electron occupies the low-energy aig 
orbital, a double degeneracy remains for the other elec- 
tron. After introducing a pseudospin-1/2 to describe the 
doubly degenerate atomic ground state, we find that their 
effective interaction is governed by a highly anisotropic 
quantum 120° Hamiltonian i^^i^^ By treating quantum 
fiuctuations using the semiclassical framework, the clas- 
sical ground state of the 120° model is shown to be stable 
against quantum fluctuations. Orbital ordering and lat- 
tice distortion derived from the classical ground states are 
consistent with the experiments. We also shown that an 
alternatively rotated orbital basis due to the staggered 
trigonal axes explains the periodic density-modulation 
observed in ab initio calculations. Moreover, since the 
very presence of trigonal distortion breaks the (approx- 
imate) rotational symmetry of t2g orbitals, orthogonal 
magnetic structure thus comes naturally from spin-orbit 
interaction and the staggering of trigonal axes. 

The rest of the paper is organized as follows. Sec. |TT] 
discusses the effective orbital 120° model on pyrochlore 
lattice and its semiclassical ground states. The corre- 
sponding orbital order and lattice distortion are discussed 
in Sec. IIIII The modulation of electron density function 
along orbital chain is addressed in Sec. IIVI The details 
of magnetic structure is presented in Sec. [V] And finally 
Sec. |VT] presents a conclusion. 



II. 120° MODEL 

The site symmetry of vanadium ions in most vana- 
dates is dominated by a cubic crystal field. Nonethe- 
less, splitting of triplet due to an additional trigonal 
distortion is known to play an important role is some 
cases. Most notably, stabilization of the unusual mag- 
netic structure in the insulating phase of V2O3 can only 
be understood when the trigonal splitting is properly 
taken into account<21"— The effects of trigonal distor- 
tions in vanadium spinels vary from one compound to 
another. For example, the trigonal splitting of t2g levels 
is essential to the understanding of heavy fermion be- 
havior in metallic LiV204.'^° On the other hand, it seems 
to have a negligible effect in another well-studied spinel 
ZnV204-^ 

Recently, experiments^ and ab initio calculation'^'^ both 
indicate a strong trigonal distortion at the vanadium 
sites of MnV204. In order to understand its effects at 
least qualitatively and to make analytical calculations 
tractable, we consider the limit of a dominating trigo- 
nal crystal field in this paper. As discussed in the in- 
troduction, the trigonal distortion still leaves a doubly 
degenerate atomic ground state. The possible long-range 
order of these localized doublets is investigated using the 
effective Hamiltonian approach. By studying the ground 




FIG. 1: Pyrochlore lattice. The numbers 0-3 denotes the 
four sublattices of pyrochlore lattice; the arrows indicate the 
local C3 axis of the respective sublattices. Explicitly, they 
are i>o = [HI], i>i = [111], i>2 = [111], and O3 = [111]. In 
120°-model Eq. ®, orbital interactions on red, green, and 
blue bonds are characterized by vectors Ayz, n^x, and iixy, 
respectively. 

state of the effective model, we discuss its implications 
for orbital and magnetic ordering in MnV204. 

In the presence of a trigonal distortion, the crystal field 
of reduced site symmetry (from cubic Oh to D^d) splits 
t2g orbitals into a singlet and a doublet separated by an 
energy gap A. The C3 symmetry axis of Z?3rf group is 
parallel to the local (111) direction of the ion (Fig. [T|). 
The aig singlet is the symmetric linear combination of 
t2g orbitals under C3 rotation 

|aig) = i^x\yz) + Vy\zx) -\- Vy\xy), (1) 

where i> = {vx-, i^yji^z) is a unit vector parallel to the local 
trigonal axis (Fig. [1]). We use the following chiral basis 
for the Cg doublet 

|e+) = i^^e-'^lyz) + i^ye+"^\zx) + i^,\xy), 

\e-) = i^^e+'^'lyz) + Vye'"^\zx) + v,\xy), ^ > 

where a; = 27r/3. A complete basis for V'^^ ion with Z(P 
configuration is given by |aige+), jaigCg), and [e^e^). 
Here the two-electron state is defined as the antisym- 
metric sum of individual one-electron states, i.e. ja/?) = 
(|a)|/3) — 1/3) |a))/'\/2. Since the aig singlet has the lowest 
energy, the atomic ground state is doubly degenerate. To 
describe the low-energy doublet manifold, we introduce 
a pseudospin-i operator r such that jr^ = ±1) are iden- 
tified with jaigC^), respectively. 

In vanadium spinels, the superexchange (SE) interac- 
tion with a 90° angle between vanadium-oxygen bonds 
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FIG. 2: Unit circle in the {T^,Ty) plane. Pseudospins par- 
ticipate in the 120° interaction (|6]) only through their projec- 
tions onto the three unit vectors Ayz, n^x, and iixv 



is dominated by direct exchange which involves electron 
hopping of the dda type^'^ 



i?SE = -5]{j2(l-SrSj)P„,,P„,j (3) 

+ (J0 S,-S, + Ji) [Pa,^{l - Pa,j) + (1 - Pa,^)Pa,J]}■ 



The various exchange constants are Jq = Jrj/{1 — 3ry), 
Ji = J(l - 77)/(l - 37?), and = J(l + + ^v), 

where J = iddal^ ^^^^ overall energy scale and rj — 
Jh/U w 0.11 denotes the ratio of Hund's exchange to 
on-site Coulomb repulsion. The subscript a = a{ij) of 
the projection operators specifies the type of orbitals in 
which electron hopping is possible between sites i and j, 
e.g. oi{ij) = xy for nearest-neighbor bonds on (110) and 
(110) chains (blue bonds in Fig. [T|). 

An interesting feature of Hamiltonian ^ is the static 
Potts-like orbital interactions which depend only on or- 
bital projection operators P^y, Pyz, and Pzx- Restricted 
to the doublet manifold jr^ = ±1), they become 
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where the three unit vectors are (Fig.[2|: 



(4) 



1 „ V3 , 



1 . 

T^zx = 2^ 



V3, 



hxy = -X. (5) 



A remark is now in order: the restricted Hilbert space 
|±) already precludes descriptions of, e.g. the ^-type 
antiferro-orbital order consisting of alternating \xy^yz) 
and \xy, zx) states;^ A more general approach is to in- 
troduce a pseudospin-1 formulation similar to the one 
used in Ref. i2Si; the trigonal splitting is then modeled by 
a spin anisotropy term. However, the resulting effective 
Hamiltonian is quite complicated and analytical calcu- 
lations are difficult. In this paper, we choose the sim- 
plified pseudospin-i formulation to explore the essential 
features of a large trigonal distortion and use perturba- 
tion method to examine the effect of excited state \e^e~) . 

The effective Hamiltonian iJ^ff of pseudospins t can be 
obtained by projecting SE Hamiltonian onto the doublet 



manifold, or equivalently by substituting the projection 
operators (HJ into Eq. ^ . In order to understand the in- 
trinsic properties of orbital interaction, we first examine 
HcS above the structural transition Ts, where nearest- 
neighbor spin correlations (Si • Sj) are isotropic. For ex- 
ample, below the Curie- Weiss temperature, the magnet is 
in a strongly correlated liquid-like state^^ the constraint 
of zero total spin Sg] = on every tetrahedron gives 



a spin correlation (S; • S, 



'S /3. In the collinear 



ferrimagnetic phase, a partial ferromagnetic order M an- 
tiparallel to the Mn moments is induced by the antiferro- 
magnetic Mn-V exchange; the spin correlation becomes 
{Si ■ Sj) = + In all cases, the effective orbital 

interaction has the following anisotropic form: 



eft 



na)(T,- • ha), 



(6) 



where Jr = |[2Ji - J2 + (S, • Sj)(J2 + 2 Jo)] > and 
a = a{ij) = yz, zx, and xy depending on the orienta- 
tion of the nearest-neighbor bond (Fig. [1} . The effective 
Hamiltonian has a form of the so-called 120° model^i2& 
which was first introduced as an effective model for per- 
ovskite eg orbital systemsi^"— Recently, the same model 
was found to describe the insulating phase of p-band 
fermions in optical lattices i^^i^^ The 120° model is closely 
related to the well-known quantum compass model^S A 
common feature shared by these highly anisotropic spin 
models is the competition between bonds along different 
directions. For compass and 120° models on bipartite 
lattices, a macroscopic degeneracy of the ground state 
results from the discrete gauge-like sliding symmetriesi^ 
Remarkably, as discussed in more detail below, such ex- 
tensive degeneracy is absent in the classical 120° model 
on non-bipartite pyrochlore lattice. 

We first discuss the origin of the macroscopic degener- 
acy in the cubic lattice 120° model. The three unit vec- 
tors fiQ. in Eq. ^ are associated with nearest-neighbor 
bonds along x, y, and z directions, respectively.^^ The 
bipartite nature of the cubic lattice allows us to trans- 
form an antiferromagnetic coupling to a ferromagnetic 
one through a 7r-rotation about axis on one sublattice. 
An unusual property of the classical model is the ap- 
pearance of planar gauge-like symmetries in addition to 
global spin rotationsi ^^i^^ A huge ground-state degener- 
acy thus results from the gauge-like Z2 transformations. 
As demonstrated in Ref. \2Q, starting from a state of uni- 
form spins, which is a ground state of the ferromagnetic 
model, another inequivalent ground state can be obtained 
by rotating all spins on a randomly chosen xy plane by 
an angle tt about Tx axis. Nonetheless, long-range or- 
der arises via the order-by-disorder mechanism which in 
general favors collinear (uniform) spin configurations.^— 

Interestingly, the above-mentioned gauge-like symme- 
try is absent for 120° model on pyrochlore lattice. This 
is because a prerequisite for the Z2 gauge-like transfor- 
mation is the existence of a subset C of lattice sites (e.g., 
planes or chains) such that pseudospins belonging to the 
subset are connected to each other by, say, either yz or 



zx bonds, while interaction of pseudospin i G C with its 
neighbor j ^ C is of the xy type exclusively. It could 
be easily checked that such a subset cannot be found 
in pyrochlore lattice. However, orbital interactions are 
still frustrated simply due to geometry: antiferromag- 
netic pseudospin interaction cannot be satisfied on all 
nearest-neighbor bonds simultaneously!^ 

Despite being geometrically frustrated, the strong 
anisotropy of 120° interaction significantly reduces the 
number of degenerate ground states. To see this, we 
note that the energy of a single bond is minimized clas- 
sically by a pair of pseodospins pointing toward ±nQ, 
respectively, where ria is the unit vector characterizing 
the anisotropic interaction of the bond. However, such 
absolute minimum can not be attained at every nearest- 
neighbor bonds due to geometrical frustration. Even 
worse, pseudospin correlation on some bonds is frus- 
trated, i.e. Ti ■ Tj > 0. In order to minimize the energy 
cost, frustrated pseudospins thus tend to align themselves 
perpendicular to ha- Through both analytical calcula- 
tion and Monte-Carlo simulations, we find that coUinear 
states with pseudospins perpendicular to either one of 
the three are the classical ground states (Fig. [3). The 
total degeneracy is six due to an additional C2-rotation 
about axis. This is in stark contrast to the macroscopic 
ground-state degeneracy of classical Heisenberg spins on 
the pyrochlore lattice.^ 

The three inequivalent ground states shown in Fig. [3] 
are characterized by the locations of frustrated bonds. 
More specifically, we introduce three staggered order 
parameters^ 

lyz = (to + Ti - T2 - T3)/4, 

hx - (t-o - Ti + r2 - T3)/4, (7) 

W = (to - Tl - r2 + T;3)/4 

to describe the orbital order. Here Ti denotes pseudospin 
average on ith sublattice. These order parameters mea- 
sure the difference of orbital configuration on bonds of 
the same type. For example, a nonzero Ixy indicates 
an antiferro-orbital order across the two xy bonds of a 
tetrahedron [see Fig. [3l^c)]. The nonzero order parame- 
ter characterizing the coUinear states of Fig. [3] are: (a) 
lyz = -^x- iy, (b) hx = ^y, and (c) l^y = y. 

This should be contrasted with the continuously degen- 
erate coUinear ground states in Heisenberg model, i.e. 
Iq, = e with e being an arbitrary unit vector;^ 

The stability of classical ground states in the pres- 
ence of quantum fiuctuations is investigated using the 
semiclassical Holstein-Primakoff transformation. We 
find that the anisotropy of the orbital exchange leads 
to a gapped quasiparticle spectrum in the whole Bril- 
louin zone. At the harmonic level, quantum fiuctua- 
tions around coUinear ground states are shown to give 
a negligible correction to the sublattice 'magnetization' 
(Ti) (about 4%), indicating the stability of the classical 
ground states. A detailed account of the semiclassical 
calculation is presented in Appendix [B] 
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(a) (b) (c) 

FIG. 3: Ground states of the orbital 120°-model. The other 
three ground states are related to the above ones by a C2 
rotation of pseudospins about the z axis. Orbital interactions 
on red, green, and blue bonds are characterized by vector Uy^, 
Uzx, and n^cy, respectively. The frustrated bonds (parallel 
pseudospins) are indicated by dashed lines. Nonzero order 
parameter characterizing the ground states are: (a) ly^ = 
-^x - iy, (b) U = -^x + |y, and (c) = y. 

Before closing this section, we remark that phonon- 
mediated orbital exchange in spinels also has the form of 
120° interaction with an effective exchange Jt oc /kp^^, 
where g is a Jahn- Teller coupling constant and kpig is the 
elastic constant of Fig phononsi^ Detailed derivation is 
presented in Appendix]^ In fact, noting that {Tx,Ty) 
forms a doublet irreducible representation of D^d group, 
the 120° type interaction in Eq. ^ is the only anisotropic 
pseudospin interaction allowed by lattice symmetry. In 
this perspective, we shall regard Jt as an effective model 
parameter in the following discussion. 

III. ORBITAL ORDER AND LATTICE 
DISTORTION 

We now discuss the orbital order and lattice distortions 
corresponding to the semiclassical ground states. Specif- 
ically, we shall focus on the coUinear state characterized 
by Ixy = +y [Fig.|Si;c)]. Since Pyz+Pzx+Pxy = 2, orbital 
orders are essentially described by linear combinations 

— P +P —IP — — T 
1 yz ' zx xy 'xj /q\ 

P2 ^ VXPzx - Pyz) ^ -Ty, 

which transform as a doublet irreducible representation 
under symmetry group D^d- Since pseudospins pointing 
along zLy directions are sitting on [110] and [110] chains, 
respectively, the ground state shown in Fig. ^c) consists 
of two distinct orbital chains characterized by Pi = 0, 
P2 = =f1i respectively. The staggered part P2 of the 
orbital order comes from the occupation difference be- 
tween yz and zx orbitals. The uniform part given by 
Pyz = Pzx = Pxy = 2/3 indicates that the three orbitals 
are equally occupied on average. 

Due to Jahn- Teller effect, a long-range orbital order 
also implies a lattice distortion in the ground state, which 
is indeed observed in MnV204 below T5 = 53 K. On 
symmetry ground, the coupling between orbital doublet 
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\tz = ±1) and distortions of the surrounding VOg octa- 
hedron has the form: 

VjT = -9iSlT^+S2Ty), (9) 

where {61,62) are coordinates of normal modes trans- 
forming as an Cg representation of group D^d- The two 
symmetry-breaking modes can be thought of as analo- 
gous to the tetragonal and orthorhombic distortions in 
a cubic VOe octahedron. The energy cost associated 
with the distortion is ^{6f + (5|), where k is an effec- 
tive elastic constant. For orbital order characterized by 
^xy = +yj minimization with respect to phonons yields 
distortions described by (5i = and 62 — ig/k on the two 
inequivalent orbital chains: octahedra on [110] and [110] 
chains are elongated along the x and y axes, respectively. 
The overall distortion preserves the tetragonal symmetry 
(lattice constants a = b > c) and is consistent with the 
observed space group lAi/ar^ 

It is interesting to note that the staggered orthorhom- 
bic distortions of VOg octahedra actually correspond to 
a softened q = lattice phonons with Fig symmetry<^i^ 
This is consistent with the fact that, by integrating 
out Fig phonons, orbital Jahn- Teller coupling gives the 
same 120° pseudospin interaction (Appendix In 
this respect, orbital ordering and structural transition in 
MnV204 can also be viewed as softening of Fig phonons 
due to cooperative Jahn- Teller effect. 

Despite the similarities between the antiferro-orbital 
order of 120° model and the A- type order proposed in 
Ref. [3, inclusion of spin-orbit interaction illustrates an 
important difference between the two cases. As we shall 
discuss later, spin-orbit coupling gives rise to an orthog- 
onal magnetic order in 120° model, which is in stark con- 
trast to collinear spins in the case of A-type orbital order. 



IV. MODULATION OF ELECTRON DENSITY 
FUNCTION 

The staggering of trigonal axes along orbital chains 
also results in a periodic variation of electron density 
distributions, despite the orbital occupation numbers are 
invariant within the chain. This is because the actual 
orbital wavefunction corresponding to t = ±y also de- 
pends on the local C3 axis: 

\Ty = +1) = V2 {v, COS ^ \X) +iyy sin C\Y)) +iy,\Z), 
\Ty = -1) = sin ^\X) + j^y cos ^\Y)) + iy,\Z). ^ 

Here the angle ^ is defined by tan ^ ( 1 — \/3) / ( 1 + \/3) , 
and \X) = \zx,xy), \Y) = \xy,yz) and \Z) = \yz,zx) are 
the two-electron basis introduced in Ref. [l3. Note that 
since Tz is diagonal in the chiral basis |±), eigenstates of 
Ty are composed of real orbitals. 

We now consider orbital chains running along [110] di- 
rection, in which the local C3 axis alternates between Uq 
and [Fig. [3]Jc)]. Along the chain, the electrons are in 




FIG. 4: Electron density of state \Ty — +1) with trigonal 
axis along (a) j>o ~ [HI] and (b) 1^3 — [111] directions. The 
explicit forms of the density corresponding to (a) and (b) are 
po(r) ± ^^(r), respectively. 

the \Ty — +1) state whose electron density can be readily 
computed 

'"W ^ IYI /^(r~r,)|(ri,r2|Tj, = +l)|'dVidV2 
= Po{r)±6p{r), (11) 

where the -I- and — signs refer to sites with Oq and trig- 
onal axes, respectively. Introducing basis functions, e.g. 
"0x^(1") = f{r)xy, where /(r) is a spherically symmet- 
ric function, the uniform and staggered parts of electron 
density are given by 

1 2 

Po(r) = -iplyir)- -cos$sin^Vy^(r)V'z^(r) (12) 
+ (2 - ^'^("^ + (2 - r^"'^) 

6p{r) = {sin^ipy^{r) + cos£,tp^a:{r)y (13) 

The resulting density functions are plotted in Fig. A 
similar density modulation, in which orbitals rotate alter- 
natively by about 45° along the orbital chain, is also ob- 
served in first-principle density functional calculations.— 
The periodic modulation of the electron density func- 
tions along orbital chains is a natural consequence of 
the staggering of local symmetry axes. In contrast, 
such density-modulation is absent in other orbital or- 
ders proposed for vanadium spinels. In A-type order, 
the two orbital chains are characterized by occupied 
two-electron states |A") and respectively. Along a 
given chain, e.g. the \X) chain, the electron density 
p{r) ~ |V'za;('') + i''/'^i/('') is invariant. In the ferro-orbital 
order proposed as the ground state of ZnV204ri there is 
only one type of orbital chain, along which the two elec- 
trons occupy states ±i\Y)) alternatively, giving 
rise to a staggered orbital angular momentum L = ±z. 
Despite the 7r-phase modulation, the electron density 

Pi^) = k^lyi^) + li'^lzi^) + V-Llr)) is the same at all 
sites. 
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FIG. 5: Magnetic orders of vanadium spinels. The and (8) 
symbols denote ±Sz components, respectively. In the limit 
A ^ A S> Jt, the spin anisotropy is dominated by Eq. (|15|) . 
and the in-plane spins point along diagonal directions Sx = 
±5y. The corresponding magnetic orders are shown in (a) 
and (b) for the case of non-magnetic and magnetic ^-site 
ions, respectively. In the opposite limit Jt 3> A 3> A, spin 
anisotropy is governed by Eq. (|17p . the resulting magnetic 
orders are shown in (c) and (d) for the respective case of non- 
magnetic and magnetic A-site ions. 



V. MAGNETIC ORDER 

In the absence of orbital order, interaction between 
vanadium spins is governed by an isotropic Heisenberg 
model, which is know to exhibit strong geometrical frus- 
tration on pyrochlore lattice>i The energy minimum of 
the model is attained by a macroscopically large number 
of states in which the total spin of every tetrahedron is 
zero = 0. The magnetic frustration is partially re- 
lieved below Tp « 56 K as the antiferromagnetic Mn-V 
exchange induces a ferrimagnetic order with antialigned 
Mn and V spins pointing along the crystal c axis. The 
transverse components of V spins remain disordered. 

The residual frustration is relieved by anisotropic spin 
exchange as well as single-ion anisotropy in the pres- 
ence of long-range orbital order. The anisotropic spin 
exchange comes from the dependence of magnetic inter- 
action on the underlying orbital configurations, as indi- 
cated by SE Hamiltonian ([3]). For example, orbital order 
corresponding to Ixy — ±y gives rise to an anisotropic ex- 
change constant such that J[iio] ~ |>^2 for bonds along 
[110] and [110] directions, i.e. directions of orbital chains, 
and J[oii] = ^[010] ~ 5|>^2 for bonds along other direc- 
tions. Consequently, upon decreasing the temperature, 
long-range antiferromagnetic spin correlation first devel- 
ops along orbital chains. However, three-dimensional 
magnetic order is not realized due to frustrated inter- 
chain couplings. 

On the other hand, single-ion anisotropy resulting from 
the relativistic spin-orbit (SO) interaction V^s — A(L-S) 
is more efficient in relieving the magnetic frustration. It 



is interesting to note that the chiral basis |e^) introduced 
in Eq. ([2]) are simultaneous eigenstates of angular mo- 
mentum operator projected onto the local trigonal axis: 



±|e 



(14) 



Restricted to the doublet subspace |±), the angular mo- 
ment operator is given by Jj = C> t^, and the effective SO 
interaction becomes 



Vls = A(S • «>) 



(15) 



We first consider limit A ^ A ^ J^. , where the single- ion 
physics dominates the Hamiltonian. The atomic ground 
state is a non-Kramers doublet (A > in V''+ ion) 



It) 



It, 



|S-i> 
|S-i> 



(16) 



Long-range ordering of spins and orbitals depends further 
on the relative strength of antiferromagnetic Mn-V and 
V-V exchanges. 

In the case of a non-magnetic A-site ion, the ground 
state consists of a uniform occupation of either | f) or 
I D states. It is easy to check that the resulting non- 
coUinear spin configuration (S.^ = +i> or = — />) also 
minimizes the V-V exchange. The corresponding ferro- 
orbital order is characterized by order parameter uIt = 
I (tq + xi + 7-2 + T3) =Fz. On the other hand, the Mn-V 
exchange is minimized by a staggering of occupied | f) 
and I D states such that the Sz component at every site 
is opposite to Mn spins (whose effect can be thought of as 
an external magnetic field). The ground state can then 
be viewed as a collection of two inequivalent spin-orbital 
chains running along [110] and [110] directions [Fig.[5jb)]. 
The corresponding orbital configuration is described by 
order parameter l^y = ±z. 

When Jt is comparable or larger than A, the competi- 
tion between SO interaction and orbital-exchange gives 
rise to a T2 cx A/J^.. The perturbation thus is essen- 
tially of order A^. To be consistent, we should take into 
account the second-order perturbations of SO interaction 
simultaneously. A straightforward calculation yields 



A^ 

VLs = ^ 



(S^ + Sy 25*^) 4- 



A2 



V^A 



Ty{Sl-Sl). (17) 



This expression can be thought of as an invariant product 
of two irreducible eg representations of D^d- Now con- 
sider ground state characterized by Ixy — +y in which 
pseudospins = ±y along [110] and [110] orbital chains, 
respectively. The second term in Eq. ([T7]) thus introduces 
a staggered spin anisotropy with easy axis parallel to y 
or X axes depending on Tj, = -f 1 or — 1, respectively. 

To determine the equilibrium spin configuration, we 
note that the two competing anisotropics Eqs. and 
(flT)) have a magnitude of order /Jt and A^/A, respec- 
tively. Consequently, when orbital exchange (including 
phonon-mediated exchange) dominates the trigonal split- 
ting Jr ^ A, the anisotropy ±{Sx — Sy) wins and aligns 
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FIG. 6: Evolution of orthogonal magnetic order as a function 
of the ratio Jt/X- The parameter and (p measure the ro- 
tation angle of pseudospin and in-plane orthogonal structure, 
respectively. The SO coupling constant is set to A = 0.75A, 
and the ferromagnetic spin component is Sz — l/\/3- The 
limiting cases (p — tv/A and (p — correspond to magnetic 
orders shown in Figs. [5](b) and (d), respectively. 



spins to either x or y directions along the respective or- 
bital chains. The resulting magnetic orders are shown 
in Figs. [5]Jc) and (d) for the case of non-magnetic and 
magnetic A-site ions, respectively. In particular, the one 
shown in Fig. El^d) is consistent with the proposed mag- 
netic ground state for MnV2 04 in Ref. [13. 

In order to understand in more detail the transition 
between these two limiting cases, we performed an ex- 
plicit calculation of the magnetic structure using Eqs. 
®, dm), and (dH). Since a detailed knowledge of Mn- 
Mn and Mn-V exchanges is required in order to com- 
pute the ferromagnetically ordered component, we set 
5"^ = 1/a/3 to simplify the calculation. As Fig. [B] shows, 
below a critical J* ~ 0.18A, pseudospins r are polarized 
along z direction, while the transverse spin components 
pointing along the diagonal directions form the orthog- 
onal structure shown in Fig. ^h). Above the critical 
J* , the transverse spins rotate uniformly (the rotation is 
described by angle (f)) while maintaining the orthogonal 
structure. At the same time, pseudospins develop a finite 
antiferro-orbital order along Ty which is characterized by 
angle tp. The calculation shows that at Jr ^ lOA, the 
spin anisotropy is already dominated by Eq. (jlZp as the 
angle cp ^ and the in-plane spins essentially point along 
either x ot y axis. 



VI. CONCLUSIONS 

To summarize, we have proposed and studied a spin- 
orbital model for vanadium spinel MnV204 taking into 
account a large trigonal distortion at the vanadium sites. 
Instead of conventional t2g triplet, our starting point is 
the doubly degenerate eg eigenstates of the trigonal crys- 
tal field. By introducing a pseudospin- 1/2 for the low- 
energy doublet, we have shown that the effective orbital 
interaction resulting from both the superexchange and 



cooperative Jahn- Teller effect is described by a quantum 
120° Hamiltonian on pyrochlore lattice. From both ana- 
lytical and numerical calculations, we have found six clas- 
sical ground states with collinear pseudospins perpendic- 
ular to either one of the three unit vectors characterizing 
the anisotropic interactions. The classical ground state is 
further shown to be stable against quantum fluctuations. 

The ground-state structure obtained from our model 
is consistent with main experimental observations and 
ab initio calculation of MnV204, namely, an antiferro- 
orbital order with tetragonal lAi/a space group, a 
density-modulation along orbital chains, and an orthogo- 
nal magnetic structure. The orbital order corresponding 
to the semiclassical ground states consists of two incquiv- 
alent orbital chains running along (110) and (110) di- 
rections, similar to the so-called A-type antiferro-orbital 
order. However, the staggering of trigonal axes along 
orbital chains gives rise to a periodic variation of the 
electron density function, which is absent in the A-type 
order. Moreover, since the trigonal distortion breaks the 
approximate rotational symmetry of t2g orbitals, orthog- 
onal magnetic structure is shown to be stabilized by the 
staggering of the single-ion spin anisotropies. 

The overall orientation of the orthogonal structure ac- 
tually depends on the relative strength of effective or- 
bital exchange Jt and spin-orbit coupling A. The ex- 
perimentally proposed orthogonal structurei^ with trans- 
verse vanadium spins pointing along either x oi y axes 
is stabilized when Jt > A. Using the large-J approach 
which assumes a dominant A over superexchange energy 
scale, we have recently shown that the same orthogo- 
nal spin structure is stabilized when the predominant 
lattice distortion is of Fig symmetry. In fact, without 
any Jahn- Teller distortion, the large-J ground state has 
a collinear antiferromagnetic order similar to the one pro- 
posed in Ref. [l3- Noting that Jt also includes contribu- 
tions from Fig phonons, results from the two comple- 
mentary approaches (large trigonal field vs large-J) are 
actually consistent: stabilization of the experimentally 
proposed orthogonal structure requires a dominant Fig 
distortion as well as a smaller spin-orbit coupling. This 
conclusion is also supported by recent ab initio calcu- 
lation which shows that inclusion of spin-orbit coupling 
does not significantly change the spin-orbital order (a fi- 
nite LS coupling, however, is still required to provide 
the spin anisotropies). Detailed analysis of Jahn- Teller 
phonons from first-principles calculation might help clar- 
ify the role of Fig phonons. 
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(a) (b) 

FIG. 7: q = lattice phonons with (a) Eg and (b) Fig symme- 
tries in a spinel structure. These normal modes are dominated 
by oxygen displacements. The red and black circles indicate 
oxygen and vanadium ions, respectively. 
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r=K2/K, 

FIG. 8: Classical phase diagram of Hamiltonian Eq. (|A1|I . 
The filled circle denotes +Tz component of pseudospin. 
The ground-state orbital order depends on the ratio r = 
K2 /Ki . The boundary separating ferro-orbital from coUinear 
antiferro-orbital states is given by Vc = 8/3. 



grating out the phonons, we obtain 



Appendix A: Phonon-mediated orbital exchange 

In this Appendix, we derive the effective pseudospin 
interaction due to cooperative Jahn- Teller effect. The 
derivation presented here is based on a related work in 
Ref. lU. To explicitly take into account the cooperative 
nature of phonon-mediated orbital exchange, we consider 
coupling of orbitals to nonlocal lattice vibrations. In fact, 
a similar study based on local Einstein-like phonons gives 
an inconclusive result i2i To further simplify the discus- 
sion, we notice that all experimentally observed struc- 
tural distortion in vanadium spinels preserves the lattice 
translational symmetry. We thus restrict our analysis to 
phonons with wavevector q = 0. 

Since orbitals mainly couple to oxygen ions, we focus 
on normal modes which are dominated by oxygen dis- 
placements. The 8 oxygen ions in a primitive unit cell of 
spinel form two tetrahedra related to each other by in- 
version symmetry. In the following, we confine ourselves 
to oxygen phonons with odd parity (coupling to even- 
parity oxygen modes cancels identically) . Among the re- 
maining modes, we find that the doublet Eg and triplet 
Fig phonons are most effective in JT coupling. The Eg 
modes correspond to the tetragonal and orthorhombic 
distortions of the oxygen tetrahedra, whereas the triply 
degenerate Fig modes represent rigid rotations of oxygen 
tetrahedra about the three cubic axes (Fig. [7]).'"^ 

As discussed in Sec. Illli the orbital doublet {Tx.,Ty) of 
trigonal crystal field couples to the Eg distortion of VOe 
octahedron [Eq. To obtain the effective orbital inter- 
action mediated by the above-mentioned lattice phonons, 
we first express the coordinates {5i,S2) describing the dis- 
tortion of a local VOe octahedron in terms of coordinates 
of lattice Eg and Fig phonons. The energy costs asso- 
ciated with the two modes are characterized by effective 
elastic constants and kp-^^, respectively. After inte- 



HjT = -Ki n ■ Tj + K^^in ■ n„)(Tj • n„), (Al) 

with the following effective exchange constants 

Ki = 2g^/kE^ + g^/kp,^ , = ig^/kp,, ■ (A2) 

The effective Hamiltonian (|A1|) contains two competing 
interactions: the Ki term denotes an isotropic Heisen- 
berg exchange, whereas the K2 term represents the 
anisotropic 120° interaction introduced in Sec. [Ill A clas- 
sical phase diagram of the above model is summarized in 
Fig. |51 where a ferro-orbital order is separated from the 
coUinear antiferro-orbital ground state of the 120° model 
discussed previously. Note that the Heisenberg term has 
a ferromagnetic sign, hence favoring a ferro-orbital order- 
ing. The ferro-orbital state is doubly degenerate with all 
pscudospins pointing along either -|-z or — z directions. 
On the other hand, psedospins in the antiferro-orbital 
phase are in a coUinear up-up-down-down configuration. 
As discussed in Sec. HH there are totally six degener- 
ate antiferro-orbital states; the one shown in Fig. [5^b) 
is characterized by a nonzero order parameter \xy = +y. 
By comparing the energies, we find a phase boundary at 
re = {K2/K1), = 8/3. 



Appendix B: Semiclassical approach to 120° model 

Here we examine quantum corrections to the classi- 
cal ground state of orbital 120° model on pyrochlore lat- 
tice. Our approach is based on a semiclassical Holstein- 
Primakoff expansion around the coUinear state described 
by order parameter l^y = +y [Fig. Wic)]. To this end, 
we first generalize the 1 20° mode l to pseudospins T of 
arbitrary length |T| = .jT{T+l): 

Hi2a^ - Jt 5Z(T, • n„)(T, • n„), (Bl) 
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FIG. 9: Quasiparticle dispersions of the 120°-model. The 
energy ek is measured in unit of Jt- The various symmetry 
points in fc-space are T = (0,0,0), K' = (1, 1,0), X = (1,0,0), 
L = a, h I), U = a, i 1), and X' = (0,0, 1). 



where Jt is an effective coupling constant. We then ex- 
pand the Hamiltonian in powers of 1 /T around the clas- 
sical ground state using Holstein-Primakoff transforma- 
tion. To simplify the calculation, we rotate the pseu- 
dospins around axis such that 



Ty = ±T. 



(B2) 



where + and — signs refer to pseudospins along [110] and 
[110] chains, respectively. The classical ground state in 
terms of rotated pseudospins is simply — +Tz. We 
then apply the standard Holstein-Primakoff transforma- 
tion: 



f+ = V2T - fltaa « a 



(B3) 



where a and satisfy the canonical boson commuta- 
tion relations. Substituting Eq. (|B3|1 into 120° Hamil- 
tonian dni), we obtain i?i2o= = -dNtJxT^ + H2 + 
where Nt is the number of unit cells, H2 of order 0{T) 
is the quadratic magnon Hamiltonian, and the omitted 
terms are of higher orders in 1/T. We note that the lin- 
ear term Hi ^ 0{T^^^) vanishes identically as expected 
for expansion around a classical ground state. 

After Fourier transformation, the quadratic Hamilto- 
nian reads 



E4 



3t' 
2 



Mk 



(B4) 



where Ok = [ao(k), • • • , a3(k), aj(— k), • • • , al{—k)]^ is an 
8-component column vector, T is a 4 x 4 identity matrix, 
and 



Mk = 



T 






Cyz 


C-zx 


^Cxy 


Cyz 





Ar 

^^xy 


Czx 




Ac 

^^xy 





Cyz 


Ac 


C-zx 


Cyz 






(B5) 



For convenience, we have defined Cap = cos[{ka -f A:^)/4] 
and Cap — cos[(fcct — ^^)/4] (the lattice constant is set 
to 1). To diagonalize H2, we consider equation of motion 
for boson operator d]^: 



.aok _ f f Z + M, 



dt 



k Mk 

-Mk 3r 



-fl-Afk 



dk = Ckdk. (B6) 



From eigenvectors of matrix £k, one can construct a 
canonical transformation Ok = 7k Ck, where matrix 7k 
satisfies 



T^r]Ti^ = r], 7^ ^Cu% = r/Ak 



(B7) 



Here the diagonal matrices rj — diag(I, —I) and 
Ak = diag(eo(k), • • • , £3(k), eo(-k), • • • ,£3(-k)). Using 
Eq. (|B7|) . the diagonalized Hamiltonian becomes 



^2 = ^ E '^kAkCk = E E ^m(k)cL(k)c™(k). (B8) 

k m k 

By setting T to its physical value 1/2, the numerically 
obtained dispersion em(k) along various symmetry di- 
rections of the Brillouin zone is shown in Fig. O The 
spectrum is fully gapped with an energy gap of order Jt, 
implying a small quantum correction. To confirm this, we 
also compute the harmonic correction to the sublattice 
'magnetization' (T^) = 1/2 — (ajai). Using the explicit 
expression of the canonical transformation 7k: 



as(k) 



,(k)c„(k)+t;,,„,(k)4(-k), (B9) 



the average quasiparticle number can be computed: 
(4 a,) = (1/Ar,)^^(at(k)a,(k)) = E™Ek<™- The 
numerically obtained sublattice 'magnetization' is given 
by (fz) w 1/2 - 0.02154. The collinear ground states 
shown in Fig. [3] thus are stable against quantum fluctu- 
ations, at least at the harmonic order. 
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